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The Linear Instability of Astrophysical Flames in Magnetic Fields 

L. J. Dursi^ 

ABSTRACT 

cn 
o . 

^ Supernovae of Type la are used as standard candles for cosmological observations despite the 

04 ' as yet incomplete understanding of their explosion mechanism. In one model, these events are 

O I thought to result from subsonic burning in the core of an accreting Carbon/Oxygen white dwarf 

that is accelerated through flame wrinkling and flame instabilities. Many such white dwarfs 
have significant magnetic fields. Here we derive the linear effects of such magnetic fields on one 
\l , flame instability, the well-known Landau-Darrieus instability. When the magnetic field is strong 

enough that the flame is everywhere sub-Alfvenic, the instability can be greatly suppressed. 
^ ■ Super- Alfvenic flames are much less affected by the field, with flames propagating parallel to 

, the fleld somewhat destabilized, and flames propagating perpendicular to the fleld somewhat 

stabilized. Trans-Alfvcnic parallel flames, however, like trans-Alfvenic parallel shocks, are seen 
■ to be non-evolutionary; understanding the behavior of these flames will require careful numerical 

simulation. 

cn 
o 

Subject headings: supernovae: general — white dwarfs - hydrodynamics — MHD — nuclear 
reactions, nucleosynthesis, abundances — conduction 



o 

to ! 1- INTRODUCTION 



X 



1.1. Supernovae Type la Progenitors 



The current standard model for Type la supernova explosions is that a flame or flames begin inside a 
' Carbon-Oxygen Chandrasekhar mass white dwarf (see for instance the review by Hillebrandt and Niemeyer 

(2000)). The burning may accelerate by flame instabilities; the Carbon-Oxygen flame is unstable to the 
Landau-Darrieus and Rayleigh- Taylor instabilities (Bychkov and Liberman 1995) which wrinkles the flame 
front, burning a larger surface area of material per unit time. 

Although the class of white dwarfs that forms the progenitors for Type la supernovae is unclear, it is 
known that many white dwarfs have signiflcant surface magnetic fields, up to 10^ G for accreting systems 
and isolated systems with magnetic fields up to 10^ G are known; central magnetic fields may be much larger 
(see the review by Wickramasinghe and Ferrario 2000, and new data in for instance Schmidt et al. 2003). 
Fields of this strength may certainly affect the dynamics of flame instabilities, and thus the mechanism of 
supernovae. It is known that magnetic flelds can destabilize shocks (see for example Stone and Edelman 1995 
and references therein), and suppress other interfacial instabilities, such as Richtmeyer-Meshkov (Samtaney 
2003) and Rayleigh-Taylor (Chandrasekhar 1981). 
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In terrestrial chemical combustion, many intermediate species are ionized, so that magnetic fields might 
be thought to have a dynamic effect. However, research on flames in magnetic fields has been sparse in the 
terrestrial flame literature because the diffusion of magnetic fields under in terrestrial burning conditions 
is fast enough that any interaction between the flame and the magnetic field is quite limited. Some recent 
work has been done on the effect of ambient magnetic fields on chemical flames (Mizutani et al. 2001; 
Wakayama and Sugie 1996) where small effects were found mainly due to the effect on large-scale motions 
of paramagnetic gases such as 02- 

In the astrophysical case magnetic fields may be more relevant, and the flame-field interactions may 
be richer in this extremely large magnetic Reynolds number regime, as perturbations in the magnetic fleld 
cannot readily diffuse away. Although there has been some look at large-scale burning dynamics in the 
presence of a large magnetic field (Ghezzi et al. 2001), little work has been presented on the microphysics of 
flame dynamics in the presence of a magnetic field. 

In this paper, we present a linear-stability analysis of the hydrodynamic flame instability, or Landau- 
Darrieus instability (Landau 1944; Darrieus 1938) in the presence of a magnetic field. 



The physical scenario of the Landau-Darrieus instability is sketched in Fig.l, where burning proceeds 
in a subsonic wave (a deflagration wave, or a flame), propagating as a plane. For simplicity, we consider 
the frame of the unperturbed flame. The flame structure itself is not modelled here; instead, we consider 
the flame to be an inflnitcsimally thin interface propagating at a known speed, S"/, into the fuel. In the 
original Landau-Darrieus derivation, the flame speed was assumed to be so much slower than the sound 
speed that the unburned fuel and burned ash regions could be assumed to have incompressible flows; in our 
astrophysical case, this assumption also holds (Timmes and Woosley 1992). We use subscript h to refer to 
quantities in the burned region, and u to refer to quantities in the unburned region. Fuel approaches the 
flame with a velocity equal to Si. 

The reaction that drives the propagating flame is assumed to be exothermic enough that the density 
behind the flame, ph is lower than the density ahead of the flame, /9„. This is expressed as a density ratio, 
= Pu/ Pb > 1- Lr this case, streamlines are refracted by the density jump across the flames, so that any 
initial wrinkling of the flame drives further wrinkling (see for instance Williams 1985, and Fig. 2.) This 
instability, driven by the density jump and mediated by hydrodynamic waves, is the well-known Landau- 
Darrieus instability, and in the linear regime it has an exponential growth rate of {e.g., Bychkov and Liberman 



where n is the growth rate, k is the wavenumber of the perturbation, and g is the gravitational acceleration. 
This instability is thought to play roles in terrestrial combustion (see for instance the references in Bychkov 
and Liberman 2000), astrophysical combustion (Hillebrandt and Niemeyer 2000), and inertial confinement 
fusion (Piriz and Portugues 2003), as well as at QCD and clcctroweak transition fronts (Kajantic 1994; Huet 
et al. 1993) with possible astrophysical application to strange stars (Benvenuto et al. 1990), polymerization 
fronts (Allah et al. 2001) and evaporation fronts such as of liquid metals (Lesin et al. 1993). 



1.2. Landau-Darrieus Instability 



2000) 




(1) 



The wrinkling perturbs the pressure and velocities on either side of the flame. If the fiame's speed into 
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the fuel is taken to be fixed at its planar value (the applicability of this assumption to astrophysical flames 
is considered in Dursi et al. 2003), the velocity at which the flame wrinkles is equal to the perturbed fluid 
velocities at the interface - e.g., the flame is advected by the perturbed velocity. 

This hydrodynamical Landau-Darrieus instability has been analytically examined in further detail in 
the astrophysical context (Bychkov and Liberman 1995; Blinnikov et al. 1995; Blinnikov and Sasorov 1996). 
In this work, we extend the derivation of the classical Landau-Darrieus instability to the case of the presence 
of a magnetic field. In §2, we consider a background magnetic field oriented along the direction of flame 
propagation. In §3, we consider a background magnetic field oriented transverse to the direction of prop- 
agation. In §4 we examine the self-consistency of assuming a fixed flame propagation speed when thermal 
diffusion will preferentially operate along magnetic field lines, In §5 we examine the effect of non-constancy 
of flame speed due to curvature, and we conclude in §6. 



2. FIELD ALONG DIRECTION OF FLAME PROPAGATION 

2.1. Perturbation equations 

Here we derive the linear growth rate of the magnetic Landau-Darrieus instability. We follow closely 
the approach and notation of Chandrasekhar (1981). We begin with the equations of incompressible mag- 
nctohydrodynamics in an inviscid, perfectly conducting fluid: 

dU^ dU., Bj fdB, dBA _ 1 dp 5p 



dt ' dxj inp \ dxj dxi ) p dxi ^ ^* p ^ ^ 



dB, _d_ 
dt dx 



{UjB,~BjU,) = (3) 



dxi 
dxi 



= (4) 
= (5) 



where Xi is the z-coordinate, Ui is the velocity in the i-coordinate direction, p is density, Sp is any fluctuation 
in the density, p is pressure, B is the magnetic field, g is the gravitational acceleration, and summation over 
repeated indicies is implied. We assume constant density states (Sp ~ 0) on the burned and unburned sides, 
with both gravity and a magnetic field pointing in the direction of propagation (which wc take here to be 
z). We consider velocity, magnetic fields, and pressure of the form 

U = {u,v,w + W) u,v,w<t:W, (6) 
B = {b^,by,h+B) b,,by,h^B, (7) 
p = P + Sp 6p-^P, (8) 

where P, B, and W are constant within each region. The velocity Wu with which the fuel approaches 
the flame is taken to be 5*;, the flame speed, so that the unperturbed flame would be motionless. We will 
assume that all velocities are highly non-relativistic, so that displacement currents and relativistic effects 
may be neglected. We also assume that the flame and flow velocities are much less than the sound speed, 
and thus may consider incompressible flow. These assumptions are all appropriate for the slow flames in 
highly degenerate material that we consider here. 
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We then expand the MHD equations to hnear order in the perturbed quantities, writing out explicitly 
the components of the vector equations. The linearization is described in more detail in Appendix A. We 
obtain: 



dt dz 

dt dz Anp dz dy ) P dy^^ ^^^^ 

(11) 



db^ 


dbA 


1 ^ R 

p dx 


dz 


dx ) 


dby 


db,\ 


1 ^ A 

pdy 


dz 


dy) 


dw 
'dt 


+ w— = 
dz 


1 ^ A 
p dz 



dbx ^^rdbx „du 



du dv dw 
dx dy dz 

dbx ^ ^ ^ ^ 
dx dy dz 



(15) 
0. (16) 



We now consider normal modes, where all perturbed quantities have a;, y, and t dependencies of the 
form exp{ikxX + ikyy + nt). In that case, and using D to denote d/dz, the equations become 

[n + WD)u^^{Dbx-ikxb,) = -'-^Sp (17) 

4:TTp P 

(n + WD)v - — (Dby -ikyb,) = -^-^Sp (18) 
inp " n p 

{n + WD)w = --D5p (19) 
P 

{n + WD)bx = BDu (20) 

{n + WD)by = BDv (21) 

{n + WD)b^ = BDw (22) 

Dw = — {ikxU + ikyv) (23) 

Db^ = - {ikxbx + ikyby) . (24) 

Adding —ikx times Eq. 17 to —iky times Eq. 18 and using Eqs. 23 and 24 gives us 

B k'^ 

(n + WD)Dw (D^-k^)b, = 5p. (25) 

inp p 

Applying D again and substituting Sp from Eq. 19 gives 

(n + WD) {D^ - k^)w - - — (D'^ - k^)Db^ = 0. (26) 

47r/9 
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Applying (n + WD) to this and applying a component of the induction equation, Eq. 22 gives 

(^{n + WDf ~ {aDf'^ {D^ - k^)w ^ (27) 

where = B'^/{4ttp) is the square of the Alfven speed. Similar equations for the other variables are: 

\n + WDf ^{aDf'jiD^ -k^)DSp = (28) 
(^{n + WDf - {aDf^ {D^ ~ k^){n + WD)b^ = 0. (29) 

We note that we are free to consider rotations around the z axis such that ky is zero and k = k^. If we do so, 
then by and v decouple from the other quantities and we need not consider them further. Then u = {i/k)Dw 
and bx = {i/k)Dbz, completing our set of equations. 

The z-depcndence of our perturbed variables, then, must be of the form 

where the C^'-* are constants within each of the burned and unburned regions, and represent the amplitudes 
of waves which travel along the z-direction with speeds — n/fc, n/fc, W , W + a, and W — a, respectively. 
Since the perturbed quantities are small (order e, as per Appendix A), these C^'^ must also be small (order 
e) quantities. 

In the case of flames propagating at the Alfven speed W ~ a in either the unburned or burned medium, 
the terms {{n + WD)^ — {aD)^) become n{n + (a + W)D), and thus only a wave corresponding to C*^^^\ 
and not C^"^ \ exists in that region. 

The unburned state will contain only waves that propagate from the flame to —oo; i.e., those with 
negative speed. Since W and a arc positive, and we are only interested in modes with the real component 
of n positive, this leaves the first wave, and the fifth if the flame is slower than the Alfven speed: Wu < au- 
The perturbations in the unburned state are then: 

5Pu = -^^^^^C(-)e^^-p„a^,Cr^e-^^ (31) 
bzu = BCi^'^e'"' +BCi'^'^e-^^^' (32) 

- + W^^fc) + a^d^-^e-^^ (33) 

k 

bxu = ^BCi^-^e^^-^B-—^ ^C(^-)e-T^^ (34) 

k{Wu - au) 

_ (n + W'^fc) au n (^-) 

" * k ^ ' k Wu-au ^ ^ ^ 

where these expressions were taken by choosing the parameters for calculating Wu from the induction 
equation and then 6p from the momentum equation, and then finding b^u from b^u using the divergence 
condition on the magnetic field and similarly using incompressibility to find ti„ from Wu- 

The same process can be repeated for the burned state. This gives us: 

dpb = p CI 'e - pba^^Cl 'e "'o+'-b - p^ai,C^ 'e »'i>-"i> (36) 

bzb = BCf^^e-''' +BCl^^^e-^^' + BCl^''>e-^^' (37) 
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b,, = -iBcf^e-^'-iB—^ d'^^'^e^w^,- _,B—^ ^d^^'e'W^^ (39) 



Note the absence of C^'^'e ("/^)^ terms; for consistency between the momentum equation and the magnetic 
field divergence condition, C^'^' must be zero. 



2.2. Jump Conditions 



Having derived the form of the perturbations in each region, we must provide the matching conditions 
at the flame, treated here as a discontinuity, to solve for the perturbations and their growth rate. Jump 
conditions in an inviscid, perfectly conducting magnetic fluid are (see for instance Shu (1992); Tidman and 
Krall (1971)): 



pUn [Ut] 

pUn [Un] -1 



Bn 

in 

V 



[Bn] 
[pUn] 

Bl 
8n 



[?7„Bt] - S„ [Ut] = 








gzf [p] 

0, 



(41) 
(42) 

(43) 

(44) 
(45) 



where [/] represents the jump in quantity / across the discontinuity, a subscript n refers to the component 
normal to the discontinuity, and subscript t refers to components tangential to the discontinuity. 

The velocity normal to the discontinuity, C/„, must be the flame speed 5*;, which is prescribed; here we 
arc considering the fixed flame speed Si = Wu ■ Since the interface is moving with respect to the far upstream 
fluid at W + Wi, where Wi is the w-velocity at the interface, the flame must be being advected — that is, 
wrinkled — with a velocity of + ui^ — [/„ = Wi . 

The directions tangential to the flame can be determined by taking the derivatives of the flame position 
Zf^ e.g.dzf / dx and dzf/dy. Because the flame position (in the unperturbed flame frame) is being wrinkled 



by the perturbed fluid velocity, and is itself a perturbed quantity Zf 
terms of Wi : 



we can express Zf m 



dt 



= — 



Wi 

n 



(46) 
(47) 



where Wi = w{zf). The difference between w{zf) and w(0) is of second order of smallness, so we consider 
Wi = w{0). This gives us the normal and two tangential directions, to linear order: 



n = 



Wi, -Wi,l 

n n 



1,0, Wi 

n 



(48) 
(49) 



We note again that we can ignore the y-component, and t2. Then we have 



Bn 

Bt 
Ut 

Un 



hx ' 

w. 



B 
ik 



-w,B 



n 
ik 



-w,W 



(51) 
(52) 

(53) 
(54) 



To zeroth order in the perturbed quantities, then, the jump conditions are: 

{B\ = 
\pW\ = 
pW [M^] + [P] = 



and hnear order gives us 



pW 



u H WiW 

n 



[Sp] 

[Wbx] - B [u] 










A further condition, 



K] = 0, 

simply states that each side of the discontinuity travels with the same velocity. 



(55) 
(56) 
(57) 
(58) 



(59) 
(60) 

(61) 
(62) 

(63) 



The continuity of bz across the interface implies (through the induction equation and the divergence-free 
property of the magnetic field) Eq. (62), so that Eq. (62) it may be discarded in favor of Eq. (59). This 
leaves us with four jump conditions: Eqns. (59), (60), (61), and (63). 

We wish to solve for the perturbed quantities, so it is the linear-order jump conditions that we are 
interested in. The zeroth-order conditions give us some restrictions on the quantities (W , pW, and B) that 
appear in the linear-order conditions. 

We apply these conditions at z = 0, as Zf itself is a perturbed quantity, and as with Wi, the difference 
between applying the boundary conditions at z = and z = Zf will be of second order in smallness. Thus, for 
any jump condition [/] = 0, we set /h(0) — /m(0) = 0. These boundary conditions can then be expressed as a 
series of linear equations in the coefficients of the z dependences, M ■ (Cu 
0. The matrix M. is 



) (^(^^) (jiA+) n(A-)\T _ 
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where n = n{kSi)~^ is the dimensionless growth rate, g ~ g{kSf)~^ is the dimensionless acceleration due to 
gravity, and a„ = auSf^ is the dimensionless Alfvcn speed in the fuel. We have expressed all quantities in 
terms of the unburned state, using p„ = apb, Si = Wu, aWu = Wb, and ^/aau = ab- 



2.3. Growth Rates 

Depending on the relation of the flame and Alfvcn velocities, different waves will be able to propagate 
away from the flame; this will determine the stability properties of the flame. The different regimes are 
shown in Fig. 3. 



2.3.1. Super-Alfvenic flames: Wu > a„; Wb > ab 

In this case, no Alfven waves can travel upstream, so that ci^ must be zero. Then we have four con- 
straints on four waves, and the system of equations represented by Eq. (64) has a solution if the determinant 
of Ai is zero; 

(Va + a„) {-a^ - + a {al + 2n)) x 
{a"^ (1 + n) + (l + n) {g -n'^) ~ a{l+g + M~2aln + gn + 3n'^ +n^)) = (65) 

Further simplification, and the removal of two trivial solutions — h ~ a ± \foLau — that lead to no actual 
wrinkling of the interface leaves us with: 

-3 _^ 1 + 3Q -2 _ + a(2a^ - g- 3) +g _ _ o? - a(g + 1) + .g ^ ^ ^^^^ 

In the limit of a„ ^ 0, this reproduces the well-known Landau-Darrieus result. 

Absent gravity, we see that the constant term in the cubic is necessarily negative, so that a positive 
real solution always exists for n — that is, the interface is always unstable. The magnetic field cannot itself 
stabilize the fiame in this case. One may examine the case of marginal stability in the presence of gravity 
by considering n = 0. This requires o? — a(g + 1) + g = 0. Since this does not include terms in Su, one 
must recover the Landau-Darrieus result; namely, that for stabilizing acceleration the condition of marginal 
stability is a = g, with larger a being unstable. 

The maximum growth rate is plotted in this case for varying g, a, and a„ in Fig. 4. We see that in this 
case, the instability is slightly enhanced over the non-magnetic (a„ = 0) case. 



2.3.2. 'Switch On' Alfvenic flames: Wu = au; Wb > ab 

When the upstream flow is exactly Alfvenic, a tangential field may "switch on" from zero in the upstream 
state to non-zero in the downstream state; such discontinuities are called "switch on" discontinuities (e.g., 
Anderson 1963). If Wu = a„, no C^f^ wave exists, as in the strictly super- Alfvenic case, and all else remains 
the same; thus, the growth rate is the a„ ^ 1 limit of the super-Alfvenic case. 
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2.3.3. Sub-Alfvenic flames: W„ < Uu; Wb < at 

In this case, upstream Alfvcn waves may propagate, and thus ' may be non-zero but C^"^ ' must be 
zero. Again, M becomes a square matrix and we may find the sohition of the hncar equations by calculating 
the determinant: 

(—1 + — n) (^a + \/aau — fi) (^Oi"^ + o?n — gh + — \fa (1 + 2a„) {g — n^) — 
a5 (1+3 + 2(2 + a„)n + fi^) + a (5(2 + 2a„ + n) + n (l + 2a^ + n2 + 2a„(l + n^ = 0. (67) 

Factoring out the two solutions (fi = a„ — 1, n = \foLa^ + a) that result in trivial solutions, one again has a 
cubic in fi, 

_3 ^1 - a + 2a„ + 2^ati _2 
1 + a 

- 2ai (2 + a„) - .g + a(l + 2a„ + 2a\ + g) _ 

H ; n 

1 + a 

(Va - l)v^ f a + + g(l + 2a„ - V^)) 

+ ^ — '- = 0. (68) 

1 + a 

Note that in this case, even absent a gravitational acceleration the constant term in the cubic is positive, 
meaning that no positive real solution need exist. 

If we express this in terms of n rather than n, and take the limit Wu — > {e.g., the interface does not 
propagate into the fluid), one can reproduce the results of the Rayleigh- Taylor instability in the presence of 
a vertical magnetic field, as discussed in § 96 of Chandrasekhar (1981). 

The scaled growth rate n as a function of a is shown in Fig. 5. The instability can be completely 
suppressed by the magnetic field. The region of stability is plotted in Fig. 6. 

2.3.4-. 'Switch Off' flames: a„ > Wu, ab = Wb 

When the downstream flow is Alfvenic, a tangential field may "switch off" in the downstream case; such 
discontinuities are called "switch off" discontinuities. 

If Wb = ab, there can be no C^^ wave, and all proceeds as in the sub-Alfvenic case with du \foL. 

2.3.5. Discussion 

Given that a sufficiently strong magnetic field suppresses the instability, it is somewhat surprising that 
while the flame is still super- Alfvenic, the growth rate of the instability actually increases with increasing 
magnetic field strength. 

The reason for this is that the presence of a magnetic field allows for a non-zero jump in the tangential 
velocity across the flame, proportional to the jump in the tangential magnetic field component. From 
the jump conditions, we can calculate the ratio of the tangential magnetic flelds across the interface as a 
function of a^; this is shown in Fig. 7. At the Alfvenic points, the flame becomes a 'switch on' or 'switch 
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off' discontinuity. For the super- Alfvcnic ffame (a„ < 1), we see ttiat tire tangential magnetic field decreases 
across the flame, and thus so docs the tangential velocity. In particular, in this case one has 

Utu^ a^af-{n + iy ^^^^ 

which is less than one. Since the magnitude of the tangential velocity decreases across the flame, the 
streamlines crossing the flame are bent still closer to the normal than in the purely hydrodynamical case 
shown in Fig. 2, and the instability is modestly enhanced. 



2.3.6. Trans- Alfvenic flames: < a^; Wi, > ai, 

In this case the coefficients C^'^ \ ct^ \ and cj^^ may all be non-zero; we have a total of three 
Alfven waves that may propagate to infinity from the flame front. Trans- Alfvenic flames then have more 
possible waves leaving the flame surface than boundary conditions to constrain them; there is no unique 
solution to the initial-value problem. 

In this case, the discontinuity is said to be non-evolutionary (see for instance Anderson 1963, or a recent 
review by Markovskii and Somov 1996.) Analytically, these discontinuities have no solution to the semi-steady 
time evolution of a small perturbation; they must then undergo some finite-sized growth or change. In the 
case of an MHD shock, the consequences of nonevolutionarity are unclear; numerical simulations of MHD 
shocks have been seen to spontaneously break up into other discontinuities (Markovskii and Skorokhodov 
2000); or stay relatively steady over long periods (Pogorelov and Matsuda 2000), depending on such things 
as shock geometry and strength. Numerical simulations will have to be performed to examine the behavior 
of flames in these circumstances. 

Simply considering the dissipative structure of an MHD shock structure does not significantly reduce the 
range of parameters in which it is non-evolutionary (Markovskii and Somov 1996). It is possible that adding 
the physics of the (compressible) flame structure to the jump conditions as was done in Pclce and Clavin 
(1982) or Matalon and Matkowsky (1982) for terrestrial non-magnetic flames may select out a particular 
solution in this case. This has never been done with a magnetic field, and is well beyond the scope of this 
work. 



It is possible, however, to constrain the possible evolution by parameterizing the allowed solutions. We 
consider taking the amplitude of the wave ci'^ ^ as a free parameter; the algebra simplifies if one considers 
a scaled ratio of this wave to the other forward-going wave. 



ia-1) 



(70) 



Note that, depending on the phase difference between the waves, 6 may be imaginary. 
In this case, one obtains the following equation for the scaled growth rate. 



_3 4a-6'_2 - a{g - 2af^ + 3) + g _ 6{g - a) 



rr + n n^ = 0, (71) 

l+a 1+a 1+a ^ ' 

which reduces (as it must) to the superalfvcnic relation if C^"^ — 0. 

One can investigate the effect of our free parameter d on the growth rate, and one finds that the 
growth rate can be greatly enhanced, if C''^"-' is very large compared to C'^"' and if they are in phase, or 
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almost stabilized if C^^"-* is very large and exactly out of phase with C^''"-'; but unless there is a significant 
stabilizing gravitational acceleration, the flame remains unstable for the entire parameter range of solutions. 
A typical case is shown in Fig. 8, with a„ = 1.3, g = 0, and a = 2, 3, 4, 5. Absent any selection criterion to 
select a unique 9, all that one can say is that the trans-alfvcnic flames are unstable, with essentially arbitrary 
growth rate. 



3. FIELD TRANSVERSE TO DIRECTION OF FLAME PROPAGATION 

3.1. Perturbation equations 

In the perpendicular-field case, the analysis is more straightforward as there is no competition between 
the flame velocity and the Alfvcn velocity. We begin with Eqs. 2-4 from §2.1, but now our magnetic field is 

bz < B. (72) 

In this case, our linearized equations in terms of normal modes are 

in + WD)u = - — Sp (73) 
P 

{n + WD)v- ^{iKby-ikyb,) = --^5p (74) 

B 1 

(n + WD)w~ - — (ik^bz-Db^) = — DSp (75) 
inp p 

{n + WD)b^ ^ ik^Bu (76) 

{n + WD)by = ik^Bv (77) 

{n + WD)b, = ik^Bw (78) 

Dw = — {ikxU + ikyv) (79) 

Dbz = — {ikxbx + ikyby) . (80) 

Going through the same steps as in § 2, but working with bx instead of b^, we obtain: 

(^{n + WDf + {ak^)^'j{D^ -k^)u = (81) 

((n + WDf + (akxf^ (D^ - k^)Sp = (82) 

(^{n + WDf + {ak^)^'j{D^ -k^){n + WD)bx = 0. (83) 

Given the linearized equations, the perturbed quantities must have z-dependencies which are linear 
combinations of the following: 

We can then find the expressions for the perturbed quantities which satisfy the momentum and induction 
equations, and the divergence criterion for the velocity and magnetic field: 

= B^C^^^^e'"' (85) 
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Wu = - - 



(n + Wuk) (l + § 



n — iabkx 



(86) 



= PJ^^cr),^^ (87) 
= _,(!l±.^i^C(-)e- (88) 



Ci'^^e'"' (89) 



= S«|^Ci^")e^^ (90) 
= -*S„^C(^")e^^ (91) 



Ufc = --^(n-M^bfc)C^'^^e-'="-a6C^^^^e"^^^^'+abC^^"^e"^^W^" (93) 
SP. - P^^^^^Cre-''^^{p,al)ci-%-"^^ (94) 



flfee "'I. (96) 



iW,B,^ e (98) 



As in § 2, C^*^' must vanish. 



3.2. Jump Conditions 

The jump conditions will be the same as in §2.2, but the components of B will be different, and both 
tangential components must be considered since one is now selected by the direction of the magnetic field. 
We have 

Bt = iB + bx,by) (99) 
B„ = —wB (100) 
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Ut = [u+—w^W,v + -^w^W 
n n 



w. 



(101) 
(102) 



Then to zeroth order in the perturbed quantities, the jump conditions give us 



and to Hnear order, 



pW 



pW [W] 



P 







[pW] = 
[WB] = 0, 



bz -WiB 

n 



u H WiW 

n 



An \ n ' 



pW 



n 



Bbr. 



Sp 



47r 







[p] 



[Wh.] = 

[Wby] = 0. 



(103) 
(104) 
(105) 

(106) 
(107) 
(108) 

(109) 

(110) 
(111) 



3.3. Growth Rates 

3.3.1. k = kyf 

Note that by Eqs. 81, 82, 83 the magnetic field only enters through ak^, so that for k^ = — that is, 
the perturbed mode is orthogonal to the magnetic field — these equations reduce to the equations for the 
Landau-Darrieus instability, and the magnetic field has no effect. 



3 . 3 . 2 . k — krf- X 

In this case, as in § 2, again the y-components decouple, and we are left with only the linear equations 
for the x-components. The linearized equations have a non-zero solution only for 

- (a^ + al- 2a^n + afi^) 
(a^ (1 + n) - al (1 + n) - a" (1 + n) (a^ + 5 + (1 + nf^ + a {2al + (1 + n) (5 - n^))) = Q. (112) 

Factoring out trivial solutions h = (a ± ia„Q;~^/^) corresponding to no wrinkling of the interface, we are left 
with 

_3 l + 3a_2 - a^{al+g + 3) + ag-al_ {a - l){a^ - a{al + g) + al) 



1 + a 



a{a + 1) 



a{l + a) 



(113) 
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In the limit of a„ ^ this reduces to the correct Landau-Darrieus hmit. Taking the other hniit. in 
unsealed quantities, of Wu 0, one obtains 

2 J pb- Pu (1 + a^) \ 

n =9k\ — — — ^ • (114) 

\Pb + Pu iirgk'^ (pb + Pu) J 

Note that the second term in this expression is different than the result quoted in Chandrasekhar (1981) § 97, 
where it is assumed that the horizontal field is equal in both regions of the domain, and so we have here 
5^(1 + a^) = + rather than 2B^. If one repeats the derivation provided therein with aBf, = i?„, as 
we have here, then the two results agree. 

The maximum real solution of Eq. (113) is plotted in Fig. 9. We see that for a„ > 1 there is significant 
suppression of growth of the instability. One can examine the stability boundary by considering n = in 
Eq. (113). This gives us 

~ {al+g)a + al^O. (115) 

The minimum a for instability as a function of a„, for varying g, is plotted in Fig. 10. Since a typical value 
for a for a carbon flame in a degenerate white dwarf is « 1.3 — 2.4 {e.g., Dursi et al. 2003), one sees that 
once Qu exceeds 1, flame instability is greatly suppressed. 

It should be noted here that the density ratio a, which in the non-magnetic case or in the case of § 2 
is given by the physicochemical properties of the burning and the EOS of the fluid, here is also in principle 
a function of the magnetic field, as the ambient field also provides pressure support. However, in the limit 
where the sound speed is much larger than the Alfven speed, this contribution is small. 



4. LOCAL EFFECTS OF MAGNETIC FIELD ON FLAME VELOCITY 

Astrophysical flames propagate by thermal diffusion. In the deep interior of a white dwarf, thermal 
diffusion is due largely to electrons; in the presence of a magnetic field, this diffusion will be anisotropic. (In 
outer regions of the star, where thermal diffusion is due to photons, there will be no such effect.) One must 
then investigate the range of validity of the assumption used in § 2 and § 3 that there will be a constant 
flame speed despite the wrinkling (and thus differing orientations) of the flame. 

We will consider here a very simple model for the anisotropic thermal diffusion, 

A..(^) = A.(0)^^±^ (116) 

where Dth is the thermal diffusivity, 9 represents the angle between the direction of diffusion — that is, 
the local direction of flame propagation — and the ambient magnetic field, and a and b parameterize the 
anisotropy. 

The local flame normal, from § 2 (but not restricted to linear order) is given by 




(117) 



Simple expressions for flame velocity {e.g., Zeldovich and Frank-Kamenetskii 1938) find that the prop- 
agation velocity scales with the square root of the diffusivity. Thus, we will consider the angle dependence 
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of the flame velocity 



4.1. Field parallel to flame propagation 

In the case of the ambient magnetic field parallel to the direction of flame propagation discussed in § 2, 

with B = Bz. We then find 

cos6i = B • n = y'l - k^zj = Vl- e^, (119) 

where Zfk = e ^ 1 as in Appendix A. 

If we then expand om' expression for Si{0) in terms of the small quantity kzf one finds 



SijO) _ b 
5,(0) 4{a + b)' 



O (e^) , (120) 



that is, there is no correction to the flame speed to the linear order of the analysis we have performed in 
§ 2. Thus, it is reasonable to consider a constant flame velocity in this case, at least for the simple angular 
dependence considered here. 



4.2. Field perpendicular to flame propagation 

On the other hand, for the case of the flame propagating across the magnetic field considered in § 3, 
B = Bx. In this case, cos 9 = 'B h = x h = —e which is of linear order, not quadratic order, in small 
quantities. In this case, we find 



^'^'^ -l + A, + o(.^) (121) 



S',(7r/2) 2a 

so that for consistency to linear order, in this case the thermal diffusion must be nearly isotropic (a 3> 6 — in 
particular, b/a must be of order e). Note that for a perpendicular field, there would have to be a significant 
isotropic component to the diffusion anyway for there to be a meaningful propagating flame, since the flame 
must propagate by diffusion across the magnetic field; however, this is a weaker condition than b/a ^ e. 
Note that for diffusivities which depend more strongly on 6, such as cos^(6'), the flame speed will no longer 
depend on linear orders of perturbed quantities even for the perpendicular field. 



5. EFFECTS OF CURVATURE 

Independent of the effects of a magnetic field, it is well known that the velocity of a laminar flame 
depends on the details of the flame structure, and that curvature of the flame will change the flame speed. 
A simple relation for the modified flame speed was first given by Markstein (1964), 

^^ = ^?(l + ^M^) (122) 

where Si is the planar flame speed. Si is the modified flame speed, and Im is an empirically determined 
'Markstein length', which is divided by the radius of curvature of the flame. This simple prescription has 
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proved very robust, and Markstcin lengths for some relevant astrophysical flames were numerically measured 
in Dursi et al. (2003). 

This non-constancy of the flame speed can be incorporated into the results of §2 and §3. One first 
observes that in the unperturbed frame, the flame now moves with a velocity 

W + w,-Si^ w, + k^ZflMW , (123) 

and thus the perturbed flame position 2:/ is now 

nzf = Wi + PzflMW (124) 

Because the flame position must be the same in both regions of the flame, the jump condition Eq. (63) 
becomes 

K] + k'^ZflM [W] = 0. (126) 

In the case of the magnetic field parallel to the direction of propagation of the flame, jump conditions 
implicitly involving Zf must be modified appropriately to use Eq. (125), and Eq. (61) has a new term directly 
from the modified flame speed. The jump conditions become 



[h] = (127) 

47r 



pW [u + ikzfW] - — [b^] = (128) 



[dp]-2k^ZflMpW[W] = gzf[p] (129) 
[Wb^]-B[u] ^ 0. (130) 

Again, the flnal jump condition may be discarded. 

One can apply these boundary conditions and flnd the solution of the linear equations as in §2. For 
super-Alfvenic flames, one obtains for the growth rate 



(1 + 3a) - 2alM _2 ("^ + a {2al - 5 - 3) + g) + 2a (a + 1) Zm _ 

^ = 0, (131) 



1 + a 1 + a 

(a - 1) (a - 5) + 2alM (c 



1 + a 

where Im = kl]\j is the non-dimensionalized Markstcin length. Note that the above expression reduces to 
Eq. (66) for hi 0. 

We can see how the Markstein length effects the growth rate in Fig. 11. If the Markstcin length is 
comparable to the wavelength of the perturbation {Mm ~ 1), the instability is greatly suppressed. 

Since the effect for other cases is similar, we don't show growth rate plots, but for completeness, we 
include the expressions for growth rate here. For sub-Alfvenic flames, using the same jump conditions one 
obtains 

_o ai - y/a (1 + 2a„) - 2aa„ + 2arA/ 

n n 

1 + a 



a 



2 



2ai(2 + a„)-.g + Q'(l + 2a„ + 2a,2 +5)_ 2aa„ (1 + V") 



n — Lm n 



1+a 1+a 
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(-1 + ^/a) y/a (^a + a2 + g ~ -Jag + 2a„g 

1 + Q' ~^ 

2a3/2 (1 + a„ - - 7^(1 +a„) + a) 



1 + a 



= 0. 



(132) 



Trans- Alfvcnic flames remain non-evolutionary, but again parameterizing the range of possible solutions 
one obtains 

^ (4q! - 6*) - 2o1m 9 

ft' -f ^ —n^ 

\ + a 

io? -f- a (a^ - g - 3) + .g) + 2a (1 + a) Im _ 
— - ^ n 

1 + a 

9{g -a) + 2alM {al - 9 - l) 
1 -(- a 



= 



(133) 



In the case of a magnetic field parallel to the flame, as in §3, the jump conditions become 

[b,-ik,ZfB] = 

pW[u + ik,zsW\~ ^{b,~ik,ZfB)[B] = 
47r 



5p - 



An 



pW[v + ikyZfW] = 

-2k^ZflMpW[W] = gzf[p\ 

[Wb,] = 

[Wby] = 0. 



(134) 
(135) 
(136) 

(137) 

(138) 
(139) 



The modifled expression for growth rate then becomes 



1 + 3a - 2q;[ 



1 + a 



M -2 

n 



-o? + iZ + al+g) -ag + al- 2lMa^ (« + 1) _ 

Tl 

a (1 + a) 

(a - 1) (a^ -a{al+ f) + al) + 2alM {a^ + al) 
a{l + a) 



(140) 



6. CONCLUSIONS 

We have examined the growth of a small perturbation to a flame in ideal MHD, and compared the 
growth to the well-known results in the absence of a field. We find that for the magnetic field to greatly 
suppress the flame instability requires a„ > Wu- This requires extremely strong magnetic fields. Assuming 
some of the largest-observed magnetic fields in white dwarfs, we have 

°" - ^ ^ ^"^ ( life) ( 10^^) (141) 

With carbon-oxygen flame speeds in the core of a white dwarf more on the order of 10® cm s~^ (Timmes and 
Woosley 1992), this is likely too low to play a major in the entire star. However, localized strong magnetic 
field regions could easily have an effect; and it is known (Cattaneo 1999) that in chaotic flows such as the 
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turbulence that is likely to exist in a rotating white dwarf that even with originally extremely weak magnetic 
fields, local peaks are quickly generated with a„ w 1. 

On the surface of a neutron star, where the magnetic field might be the same but densities would be 
lower, 



and magnetic field effects might well be significant everywhere. 

Weaker magnetic fields also effect flame growth, but only modestly. If the fiame is everywhere super- 
Alfvenic, the effects of magnetic field will somewhat stabilize the fiame if the magnetic field is perpendicular 
to the flame propagation, and destabilize it if it is parallel. 

It is instructive to compare the results presented here with those known for MHD shocks, where fast 
parallel shocks (which are super- Alfvenic) and perpendicular shocks are stable under a broad range of 
parameters (Gardner and Kruskal 1964), while slow shocks (which are sub- Alfvenic) are broadly unstable 
(Lesson and Deshpande 1967; Edel'man 1990; Stone and Edelman 1995); this is just the opposite of what is 
reported here. The key difference between the two, for the purposes of this instability, is the direction of the 
density jump across the discontinuity — recall that it is this density jump which drives the instability. 

In the flames discussed here, heating from the reactions ensures that pf, < p„; that is, a > 1. However, 
in the case of a shock, in our notation a < 1. This changes the underlying driving of the instability; indeed, 
if one uses a < 1 in the growth rates presented above, one flnds super- Alfvenic and parallel flames to be 
stable, and sub- Alfvenic flames to be unstable. 

This difference also changes the properties of the nonevolutionary nature of the interface in the trans- 
Alfvenic case. In the case of a trans-Alfvcnic MHD shock, there is one too few Alfven waves to satisfy the 
boundary conditions, so there is no solution; in the case of a trans-Alfvcnic flame, there is one too many, 
so there are an inflnite number of solutions, and no unique solution. This will change the nature of the 
behavior in this regime. Future work must include numerical simulations of the interesting trans-Alfvcnic 
flame regime, which could occur locally in a white dwarf and in large regions on the surface of a neutron 
star. 

This work was supported by the Department of Energy Computational Science Graduate Fellowship 
Program of the Office of Scientific Computing and Office of Defense Programs in the Department of Energy 
under contract DE-FG02-97ER25308. The author thanks R. Rosner, A. Calder, and T. Dupont for helpful 
comments, and R. Morgan-Dursi for her support. 
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A. 



Linearization 



A.l. MHD equations 



In this section, we describe in more detail our process of linearizing the ideal MHD equations. 



We perturb the initially planar fiame with a very small amplitude wrinkle, giving it an initial perturbed 
position Zf. We require the amplitude to be small compared to the wavelength of the perturbation: 
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— i=z/fc = e<l. (Al) 
A 

Such a perturbation will generate perturbations in the magnetic fields and velocities of the same order; 
we expand them in series of this perturbation parameter e. In the case of a parallel magnetic field, 

B = B(e6^.i+e26^2 + ---,efeyi+e^&y2 + ---'l + e^-i+e^^-2 + ---) (A2) 
U = VF {eui + e^U2 H ,evi + (?v-2 H , 1 + ewi + H ) . (A3) 

Because the flow is incompressible, the pressure fluctuations will be fluctuations in the dynamic, not 
thermodynamic, pressure; thus the pressure fluctuations scale with pW^ rather than P: 

p = P + pW^eSpi + pW^e^6p2 + ■■■■ (A4) 

We assume that all the non-dimensional expansion coefficients are of order unity. We will also assume 
that the non-dimensional parameters of the problem — the scaled growth rate (n = n/{kW)), gravity 
(5 ~ .9/(^^^)) ^-i^d Alfvcn speed = a^/Wu) — are order unity quantities. 

We may now write the ideal MHD equations in terms of these expansions. W, B, and e are independent 
of time or position, and we retain terms up to first order in e. For the momentum equation, we then have: 

Re-expressing the linear-order perturbed quantities in dimensional units, e.g. {u, v, w) — eW{ui, vi, Wi), 
{bx,by,bz) = eB{bri.i,by^,bzi), Sp = epW^Spi, we obtain Eqns. (9-11). Other equations, and the boundary 
conditions, arc linearized in the same manner. 



A. 2. Matching Boundary Conditions 

In applying the boundary conditions at z = Zf, we can consider expanding the exponential z-depcndencies 
in Taylor series. In the case of the magnetic field, expanding the exponential terms on the unburned side 
(the terms on the burned side are analogous) one obtains: 



1 - e + + • • • (A8) 



n 



2 



= I - —l^e + ——^e^ + ■ ■ ■ (A9) 
1 + au (1 + a^y 

e-w^^f = l--^e+ ^, e' + ---, (AlO) 

1 - flu [1 - aii^y 

while the values for at z = would clearly be exactly unity. Since these terms multiply the already order-e 
coefficients C*-*\ then, the difference in the perturbed quantities between their values at z = and z = Zf is 
of order and can be neglected in this analysis as long as 1 — a„ 76 e. As we will see, however, as a„ 1, 
the terms of the form (AlO) vanish, removing this restriction, and we may consider zj — for the boundary 
conditions. 
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flame 

Fig. 1. — A sketch of the physical setup of the Landau-Darrieus instability in the frame of the flame front. 
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Unburned 




Burned 



Fig. 2. — Because the tangential component of velocity is constant across the flame in the hydrodynamical 
case, while the normal component jumps, streamlines crossing the flame are bent towards the flame normal. 
Since at large distances the streamlines must remain parallel, the stream lines 'fan out' near flame surfaces 
curved towards the unburned gas; this locally reduces the flow velocity, allowing the flame in that region to 
propagate further ahead. (Figure taken after Williams (1985) Fig. 9.8). 
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super- Alfvenic 



trans- Alfvenic 



Fig. 3.^ — ' The waves and their speeds propagating from the flame (vertical line), following Anderson (1963), 
in the case of the field parallel to flame propagation. The waves shown, top to bottom, are C^^ ^ /C^^ \ 
and C^'^ From left to right, the states are: super- Alfvenic (a^ < 1); 'switch on' Alfvenic, 
{qu = 1); trans- Alfvenic (1 < < ^/a); 'switch off' Alfvenic, (a„ = y/a); and sub-Alfvenic (a„ > ^/a). 



2.D . 2.0 2.0 



1,5 ■ 1.5- 1.5 




Fig. 4. — The growth rate for a parallel magnetic field when the flame is everywhere super- Alfvenic. We 
plot the maximum real part of the scaled growth rate (n) as a function of a„, with g ~ —1 (destabilizing) 
on the left, g = in the center, and g = +1 (stabilizing) on the right. Lines are plotted for. top to bottom, 
a = 8, 4, 2, 1.5, 1.25, 1.1. 
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Fig. 5. — The growth rate for a paraUel magnetic field in the sub-Alfvcnie case with destabihzing gravity; 
with zero or stabihzing gravity, the flame is stable. We plot the maximum real part of the scaled growth 
rate (n) as a function of a„, with g = —1. Lines are plotted for, left to right, a = 1.1, 2, 4, 6. 
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Fig. 6. — Stability boundaries in a for the case of a parallel magnetic field as a function of Su. Note that 
for the superalfvcnic case, the region above the curves is unstable, and for the subalfvcnic case, the region 
below the curves is unstable. 
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a 




Fig. 7. — Plotted is the ratio of the tangential magnetic field at the interface in the burned fluid to that in 
the fuel, Bth/Btu = ([-Bt]/-Sf„) + 1. At a.,; = 1, the flame is a 'switch off' discontinuity; at a„ — ^/a, the 
flame is a 'switch on' discontinuity. 
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-20 -10 10 20 -20 -15 -10 -5 5 10 15 20 

Fig. 8. — On the left, the maximum real part of the scaled growth rate (n) in the trans- Alfvenic case, with 
g — 0, du = 1.3, and a (top to bottom) 5,4,3,2. On the right, the corresponding imaginary component. For 
any value of 9 the flame is unstable; the flame can only partially stabilize itself by sending forward fairly 
strong (several times the amplitude of the C*-'*"-' wave) and oscillatory Alfven waves. 
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Fig. 9. — The maximum real part of the scaled growth rate (n) for the case of a perpendicular magnetic 
field as a function of a^, with g = —1 on the left, g = in the center, and g = +1 on the right. Lines are 
plotted for, top to bottom, a ~ 4, 2, 1.5, 1.25, and 1.1 



a 




Fig. 10. — Stability boundaries in a for the case of a perpendicular magnetic field as a function of a„, for 
(top to bottom) g = 5,2,1,0,-1,-2. 
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Fig. 11. — The growth rate for a paraUel magnetie field when the flame is everywhere super- Alfvcnic, and 
has a negative Markstein length. We plot the maximum real part of the scaled growth rate (n) as a function 
of Qu, with Im = on the top left, Im = —0.05 on top right, Im — —0.10 on the bottom left, and Im — —0.25 
on the bottom right. Lines are plotted for, top to bottom, a = 8, 4, 2, 1.5, 1.25, 1.1. 



